/*=========================================================================
Elite College as Engines of Upward Mobility: Evidence from Colombia's Ser Pilo Paga
Authors: Juliana Londoño-Vélez, Catherine Rodriguez, Fabio Sánchez
and Luis Esteban Álvarez-Arango

Creation date: June 6th, 2025
--------------------------------------------------------------------------
Figures C1a, C1b, C2a, C2b, C3a, C3b, C4a, C4b, C5a and C5b
=========================================================================*/

********************************************************************************
*** HETEROGENEITY: ACCESS_0_HQ
********************************************************************************

use 	"${Intermediate}/heterogeneities_access_0_hq.dta", clear

replace nombre = "cc" if _n==1
replace nombre = "cr" if _n==2
replace nombre = "sec" if _n==3
replace nombre = "ser" if _n==4
replace nombre = "cs" if _n==5
replace nombre = "ci" if _n==6
replace nombre = "cc0" if _n==7
replace nombre = "cr0" if _n==8
replace nombre = "pc" if _n==9
replace nombre = "pr" if _n==10
replace nombre = "bwc" if _n==11
replace nombre = "bwr" if _n==12
replace nombre = "nec0" if _n==13
replace nombre = "nec" if _n==14
replace nombre = "ner0" if _n==15
replace nombre = "ner" if _n==16
replace nombre = "n" if _n==17

xpose, clear varname

rename v1 cc
rename v2 cr
rename v3 sec
rename v4 ser
rename v5 cs
rename v6 ci
rename v7 cc0
rename v8 cr0
rename v9 pc
rename v10 pr
rename v11 bwc
rename v12 bwr
rename v13 nec0
rename v14 nec
rename v15 ner0
rename v16 ner
rename v17 n
rename _varname nombre
drop if _n==1

order nombre
		
* 		Plotting conventional RD coefficient and conventional confidence intervals
drop 	cs ci
gen 	cs = cc+(1.96*sec)
gen 	ci = cc-(1.96*sec)

gen 	p = _n
replace p = _n-33 if p>33

gen 	r = "11" if regexm(nombre,"_saber11")
replace r = "si" if regexm(nombre,"_sisben")

replace nombre = subinstr(nombre, "_saber11", "",.)
replace nombre = subinstr(nombre, "_sisben", "",.)
replace nombre = "Baseline" if nombre=="total"
replace nombre = strproper(nombre)
replace nombre = "Not ethnic" if nombre=="Noethnia"
replace nombre = "Ethnic" if nombre=="Ethnia"
replace nombre = "First-gen college" if nombre=="Noparental_Ba"
replace nombre = "Not first-gen college" if nombre=="Parental_Ba"
replace nombre = "HS HQ college <P50" if nombre=="Enrollment_Hq_1"
replace nombre = "HS HQ college >P50" if nombre=="Enrollment_Hq_2"
replace nombre = "HS SABER 11 Q1" if nombre=="Quartil_Sb11_1"
replace nombre = "HS SABER 11 Q2" if nombre=="Quartil_Sb11_2"
replace nombre = "HS SABER 11 Q3" if nombre=="Quartil_Sb11_3"
replace nombre = "HS SABER 11 Q4" if nombre=="Quartil_Sb11_4"
replace nombre = "HS Higher Ed Q1" if nombre=="Quartil_Enrollment_1"
replace nombre = "HS Higher Ed Q2" if nombre=="Quartil_Enrollment_2"
replace nombre = "HS Higher Ed Q3" if nombre=="Quartil_Enrollment_3"
replace nombre = "HS Higher Ed Q4" if nombre=="Quartil_Enrollment_4"
replace nombre = "HS P310 Q1" if nombre=="Quartil_P310_1"
replace nombre = "HS P310 Q2" if nombre=="Quartil_P310_2"
replace nombre = "HS P310 Q3" if nombre=="Quartil_P310_3"
replace nombre = "HS P310 Q4" if nombre=="Quartil_P310_4"
replace nombre = "Not selective college" if nombre=="College_Selectivity_1"
replace nombre = "Selective college" if nombre=="College_Selectivity_2"
replace nombre = "Not selective program" if nombre=="Major_Selectivity_1"
replace nombre = "Selective program" if nombre=="Major_Selectivity_2"

drop 	if  regexm(nombre, "Selective") | regexm(nombre, "Not selective")
drop 	if  regexm(nombre, "P310")
drop 	if  regexm(nombre, "HS Higher Ed")
drop 	if  regexm(nombre, "Tercil") | regexm(nombre, "Stratum") 

replace p = 2 if regexm(nombre, "Male")
replace p = 3 if regexm(nombre, "Female")
replace p = 4 if regexm(nombre, "Not ethnic")
replace p = 5 if regexm(nombre, "Ethnic")
replace p = 6 if regexm(nombre, "First-gen college")
replace p = 7 if regexm(nombre, "Not first-gen college")
replace p = 8 if regexm(nombre, "Rural")
replace p = 9 if regexm(nombre, "Urban")
replace p = 10 if nombre == "HS SABER 11 Q1" 
replace p = 11 if nombre == "HS SABER 11 Q2"
replace p = 12 if nombre == "HS SABER 11 Q3" 
replace p = 13 if nombre == "HS SABER 11 Q4" 
replace p = 14 if nombre == "HS HQ college <P50" 
replace p = 15 if nombre == "HS HQ college >P50"

drop 	if p>15
labmask p, values(nombre)

twoway 	(rcap cs ci p if r=="11", lcolor(black)) ///
		(scatter cc p if r=="11", ysc(r(0)) mcolor(black) msymbol(O)) ///
		,ytitle(RD-IV Coefficient and 95% CI) yline(0) xlab(1(1)15,valuelabel angle(45)) ///
		legend(off) xsize(8) xtitle("") ysc(r(0 1)) ylab(0(.2)1) ///
		caption(Sample restricted to SISBEN-eligible individuals, size(small)) graphregion(color(white)) 
		graph export "${outcomes}/FigC1a.pdf", replace width(20)
	
twoway 	(rcap cs ci p if r=="si", lcolor(black)) ///
		(scatter cc p if r=="si", ysc(r(0)) mcolor(black) msymbol(O)) ///
		,ytitle(RD-IV Coefficient and 95% CI) yline(0) xlab(1(1)15,valuelabel angle(45))  ///
		legend(off) xsize(8) xtitle("") ysc(r(0 1)) ylab(0(.2)1) ///
		caption(Sample restricted to SABER 11-eligible individuals, size(small)) graphregion(color(white)) 
		graph export "${outcomes}/FigC1b.pdf", replace width(20)

********************************************************************************
*** HETEROGENEITY: GRADUATION_EXAM_PRO_HQ
********************************************************************************		
use 	"${Intermediate}/heterogeneities_graduation_exam_pro_hq.dta", clear

replace nombre = "cc" if _n==1
replace nombre = "cr" if _n==2
replace nombre = "sec" if _n==3
replace nombre = "ser" if _n==4
replace nombre = "cs" if _n==5
replace nombre = "ci" if _n==6
replace nombre = "cc0" if _n==7
replace nombre = "cr0" if _n==8
replace nombre = "pc" if _n==9
replace nombre = "pr" if _n==10
replace nombre = "bwc" if _n==11
replace nombre = "bwr" if _n==12
replace nombre = "nec0" if _n==13
replace nombre = "nec" if _n==14
replace nombre = "ner0" if _n==15
replace nombre = "ner" if _n==16
replace nombre = "n" if _n==17

xpose, clear varname

rename v1 cc
rename v2 cr
rename v3 sec
rename v4 ser
rename v5 cs
rename v6 ci
rename v7 cc0
rename v8 cr0
rename v9 pc
rename v10 pr
rename v11 bwc
rename v12 bwr
rename v13 nec0
rename v14 nec
rename v15 ner0
rename v16 ner
rename v17 n
rename _varname nombre
drop if _n==1

order nombre
            
* 		Plotting conventional RD coefficient and conventional confidence intervals
drop 	cs ci
gen 	cs = cc+(1.96*sec)
gen 	ci = cc-(1.96*sec)

gen 	p = _n
replace p = _n-33 if p>33

gen 	r = "11" if regexm(nombre,"_saber11")
replace r = "si" if regexm(nombre,"_sisben")

replace nombre = subinstr(nombre, "_saber11", "",.)
replace nombre = subinstr(nombre, "_sisben", "",.)
replace nombre = "Baseline" if nombre=="total"
replace nombre = strproper(nombre)
replace nombre = "Not ethnic" if nombre=="Noethnia"
replace nombre = "Ethnic" if nombre=="Ethnia"
replace nombre = "First-gen college" if nombre=="Noparental_Ba"
replace nombre = "Not first-gen college" if nombre=="Parental_Ba"
replace nombre = "HS HQ college <P50" if nombre=="Enrollment_Hq_1"
replace nombre = "HS HQ college >P50" if nombre=="Enrollment_Hq_2"
replace nombre = "HS SABER 11 Q1" if nombre=="Quartil_Sb11_1"
replace nombre = "HS SABER 11 Q2" if nombre=="Quartil_Sb11_2"
replace nombre = "HS SABER 11 Q3" if nombre=="Quartil_Sb11_3"
replace nombre = "HS SABER 11 Q4" if nombre=="Quartil_Sb11_4"
replace nombre = "HS Higher Ed Q1" if nombre=="Quartil_Enrollment_1"
replace nombre = "HS Higher Ed Q2" if nombre=="Quartil_Enrollment_2"
replace nombre = "HS Higher Ed Q3" if nombre=="Quartil_Enrollment_3"
replace nombre = "HS Higher Ed Q4" if nombre=="Quartil_Enrollment_4"
replace nombre = "HS P310 Q1" if nombre=="Quartil_P310_1"
replace nombre = "HS P310 Q2" if nombre=="Quartil_P310_2"
replace nombre = "HS P310 Q3" if nombre=="Quartil_P310_3"
replace nombre = "HS P310 Q4" if nombre=="Quartil_P310_4"
replace nombre = "Not selective college" if nombre=="College_Selectivity_1"
replace nombre = "Selective college" if nombre=="College_Selectivity_2"
replace nombre = "Not selective program" if nombre=="Major_Selectivity_1"
replace nombre = "Selective program" if nombre=="Major_Selectivity_2"

drop 	if  regexm(nombre, "Selective") | regexm(nombre, "Not selective")
drop 	if  regexm(nombre, "P310")
drop 	if  regexm(nombre, "HS Higher Ed")
drop 	if  regexm(nombre, "Tercil") | regexm(nombre, "Stratum") 

replace p = 2 if regexm(nombre, "Male")
replace p = 3 if regexm(nombre, "Female")
replace p = 4 if regexm(nombre, "Not ethnic")
replace p = 5 if regexm(nombre, "Ethnic")
replace p = 6 if regexm(nombre, "First-gen college")
replace p = 7 if regexm(nombre, "Not first-gen college")
replace p = 8 if regexm(nombre, "Rural")
replace p = 9 if regexm(nombre, "Urban")
replace p = 10 if nombre == "HS SABER 11 Q1" 
replace p = 11 if nombre == "HS SABER 11 Q2"
replace p = 12 if nombre == "HS SABER 11 Q3" 
replace p = 13 if nombre == "HS SABER 11 Q4" 
replace p = 14 if nombre == "HS HQ college <P50" 
replace p = 15 if nombre == "HS HQ college >P50"

drop 	if p>15
labmask p, values(nombre)

twoway 	(rcap cs ci p if r=="11", lcolor(black)) ///
		(scatter cc p if r=="11", ysc(r(0)) mcolor(black) msymbol(O)) ///
		,ytitle(RD-IV Coefficient and 95% CI) yline(0) xlab(1(1)15,valuelabel angle(45)) ///
		legend(off) xsize(8) xtitle("") ylab(0(.1).7) ///
		caption(Sample restricted to SISBEN-eligible individuals, size(small)) graphregion(color(white)) 
		graph export "${outcomes}/FigC2a.pdf", replace width(20)
	
twoway 	(rcap cs ci p if r=="si", lcolor(black)) ///
		(scatter cc p if r=="si", ysc(r(0)) mcolor(black) msymbol(O)) ///
		,ytitle(RD-IV Coefficient and 95% CI) yline(0) xlab(1(1)15,valuelabel angle(45))  ///
		legend(off) xsize(8) xtitle("") ylab(-.4(.2).8) ///
		caption(Sample restricted to SABER 11-eligible individuals, size(small)) graphregion(color(white)) 
		graph export "${outcomes}/FigC2b.pdf", replace width(20)

********************************************************************************
*** HETEROGENEITY: GRADUATION_EXAM_PRO
********************************************************************************
use 	"${Intermediate}/heterogeneities_score_saberpro_5.dta", clear

replace nombre = "cc" if _n==1
replace nombre = "cr" if _n==2
replace nombre = "sec" if _n==3
replace nombre = "ser" if _n==4
replace nombre = "cs" if _n==5
replace nombre = "ci" if _n==6
replace nombre = "cc0" if _n==7
replace nombre = "cr0" if _n==8
replace nombre = "pc" if _n==9
replace nombre = "pr" if _n==10
replace nombre = "bwc" if _n==11
replace nombre = "bwr" if _n==12
replace nombre = "nec0" if _n==13
replace nombre = "nec" if _n==14
replace nombre = "ner0" if _n==15
replace nombre = "ner" if _n==16
replace nombre = "n" if _n==17

xpose, clear varname

rename v1 cc
rename v2 cr
rename v3 sec
rename v4 ser
rename v5 cs
rename v6 ci
rename v7 cc0
rename v8 cr0
rename v9 pc
rename v10 pr
rename v11 bwc
rename v12 bwr
rename v13 nec0
rename v14 nec
rename v15 ner0
rename v16 ner
rename v17 n
rename _varname nombre
drop if _n==1

order nombre
            
* 		Plotting conventional RD coefficient and conventional confidence intervals
drop 	cs ci
gen 	cs = cc+(1.96*sec)
gen 	ci = cc-(1.96*sec)

gen 	p = _n
replace p = _n-33 if p>33

gen 	r = "11" if regexm(nombre,"_saber11")
replace r = "si" if regexm(nombre,"_sisben")

replace nombre = subinstr(nombre, "_saber11", "",.)
replace nombre = subinstr(nombre, "_sisben", "",.)
replace nombre = "Baseline" if nombre=="total"
replace nombre = strproper(nombre)
replace nombre = "Not ethnic" if nombre=="Noethnia"
replace nombre = "Ethnic" if nombre=="Ethnia"
replace nombre = "First-gen college" if nombre=="Noparental_Ba"
replace nombre = "Not first-gen college" if nombre=="Parental_Ba"
replace nombre = "HS HQ college <P50" if nombre=="Enrollment_Hq_1"
replace nombre = "HS HQ college >P50" if nombre=="Enrollment_Hq_2"
replace nombre = "HS SABER 11 Q1" if nombre=="Quartil_Sb11_1"
replace nombre = "HS SABER 11 Q2" if nombre=="Quartil_Sb11_2"
replace nombre = "HS SABER 11 Q3" if nombre=="Quartil_Sb11_3"
replace nombre = "HS SABER 11 Q4" if nombre=="Quartil_Sb11_4"
replace nombre = "HS Higher Ed Q1" if nombre=="Quartil_Enrollment_1"
replace nombre = "HS Higher Ed Q2" if nombre=="Quartil_Enrollment_2"
replace nombre = "HS Higher Ed Q3" if nombre=="Quartil_Enrollment_3"
replace nombre = "HS Higher Ed Q4" if nombre=="Quartil_Enrollment_4"
replace nombre = "HS P310 Q1" if nombre=="Quartil_P310_1"
replace nombre = "HS P310 Q2" if nombre=="Quartil_P310_2"
replace nombre = "HS P310 Q3" if nombre=="Quartil_P310_3"
replace nombre = "HS P310 Q4" if nombre=="Quartil_P310_4"
replace nombre = "Not selective college" if nombre=="College_Selectivity_1"
replace nombre = "Selective college" if nombre=="College_Selectivity_2"
replace nombre = "Not selective program" if nombre=="Major_Selectivity_1"
replace nombre = "Selective program" if nombre=="Major_Selectivity_2"

drop 	if  regexm(nombre, "Selective") | regexm(nombre, "Not selective")
drop 	if  regexm(nombre, "P310")
drop 	if  regexm(nombre, "HS Higher Ed")
drop 	if  regexm(nombre, "Tercil") | regexm(nombre, "Stratum") 

replace p = 2 if regexm(nombre, "Male")
replace p = 3 if regexm(nombre, "Female")
replace p = 4 if regexm(nombre, "Not ethnic")
replace p = 5 if regexm(nombre, "Ethnic")
replace p = 6 if regexm(nombre, "First-gen college")
replace p = 7 if regexm(nombre, "Not first-gen college")
replace p = 8 if regexm(nombre, "Rural")
replace p = 9 if regexm(nombre, "Urban")
replace p = 10 if nombre == "HS SABER 11 Q1" 
replace p = 11 if nombre == "HS SABER 11 Q2"
replace p = 12 if nombre == "HS SABER 11 Q3" 
replace p = 13 if nombre == "HS SABER 11 Q4" 
replace p = 14 if nombre == "HS HQ college <P50" 
replace p = 15 if nombre == "HS HQ college >P50"

drop 	if p>15
labmask p, values(nombre)

twoway 	(rcap cs ci p if r=="11", lcolor(black)) ///
		(scatter cc p if r=="11", ysc(r(0)) mcolor(black) msymbol(O)) ///
		,ytitle(RD-IV Coefficient and 95% CI) yline(0) xlab(1(1)15,valuelabel angle(45)) ///
		legend(off) xsize(8) xtitle("") ylab(-.2(.2).8) ///
		caption(Sample restricted to SISBEN-eligible individuals, size(small)) graphregion(color(white)) 
		graph export "${outcomes}/FigC3a.pdf", replace width(20)
	
twoway 	(rcap cs ci p if r=="si", lcolor(black)) ///
		(scatter cc p if r=="si", ysc(r(0)) mcolor(black) msymbol(O)) ///
		,ytitle(RD-IV Coefficient and 95% CI) yline(0) xlab(1(1)15,valuelabel angle(45))  ///
		legend(off) xsize(8) xtitle("") ylab(-1(.5)2) ///
		caption(Sample restricted to SABER 11-eligible individuals, size(small)) graphregion(color(white)) 
		graph export "${outcomes}/FigC3b.pdf", replace width(20)

********************************************************************************
*** HETEROGENEITY: WAGE
********************************************************************************
	
use 	"${Intermediate}/heterogeneities_mw_9.dta", clear

replace nombre = "cc" if _n==1
replace nombre = "cr" if _n==2
replace nombre = "sec" if _n==3
replace nombre = "ser" if _n==4
replace nombre = "cs" if _n==5
replace nombre = "ci" if _n==6
replace nombre = "cc0" if _n==7
replace nombre = "cr0" if _n==8
replace nombre = "pc" if _n==9
replace nombre = "pr" if _n==10
replace nombre = "bwc" if _n==11
replace nombre = "bwr" if _n==12
replace nombre = "nec0" if _n==13
replace nombre = "nec" if _n==14
replace nombre = "ner0" if _n==15
replace nombre = "ner" if _n==16
replace nombre = "n" if _n==17

xpose, clear varname

rename v1 cc
rename v2 cr
rename v3 sec
rename v4 ser
rename v5 cs
rename v6 ci
rename v7 cc0
rename v8 cr0
rename v9 pc
rename v10 pr
rename v11 bwc
rename v12 bwr
rename v13 nec0
rename v14 nec
rename v15 ner0
rename v16 ner
rename v17 n
rename _varname nombre
drop if _n==1

order nombre
            
*** 		Plotting conventional RD coefficient and conventional confidence intervals
drop 	cs ci
gen 	cs = cc+(1.96*sec)
gen 	ci = cc-(1.96*sec)

gen 	p = _n
replace p = _n-33 if p>33

gen 	r = "11" if regexm(nombre,"_saber11")
replace r = "si" if regexm(nombre,"_sisben")

replace nombre = subinstr(nombre, "_saber11", "",.)
replace nombre = subinstr(nombre, "_sisben", "",.)
replace nombre = "Baseline" if nombre=="total"
replace nombre = strproper(nombre)
replace nombre = "Not ethnic" if nombre=="Noethnia"
replace nombre = "Ethnic" if nombre=="Ethnia"
replace nombre = "First-gen college" if nombre=="Noparental_Ba"
replace nombre = "Not first-gen college" if nombre=="Parental_Ba"
replace nombre = "HS HQ college <P50" if nombre=="Enrollment_Hq_1"
replace nombre = "HS HQ college >P50" if nombre=="Enrollment_Hq_2"
replace nombre = "HS SABER 11 Q1" if nombre=="Quartil_Sb11_1"
replace nombre = "HS SABER 11 Q2" if nombre=="Quartil_Sb11_2"
replace nombre = "HS SABER 11 Q3" if nombre=="Quartil_Sb11_3"
replace nombre = "HS SABER 11 Q4" if nombre=="Quartil_Sb11_4"
replace nombre = "HS Higher Ed Q1" if nombre=="Quartil_Enrollment_1"
replace nombre = "HS Higher Ed Q2" if nombre=="Quartil_Enrollment_2"
replace nombre = "HS Higher Ed Q3" if nombre=="Quartil_Enrollment_3"
replace nombre = "HS Higher Ed Q4" if nombre=="Quartil_Enrollment_4"
replace nombre = "HS P310 Q1" if nombre=="Quartil_P310_1"
replace nombre = "HS P310 Q2" if nombre=="Quartil_P310_2"
replace nombre = "HS P310 Q3" if nombre=="Quartil_P310_3"
replace nombre = "HS P310 Q4" if nombre=="Quartil_P310_4"
replace nombre = "Not selective college" if nombre=="College_Selectivity_1"
replace nombre = "Selective college" if nombre=="College_Selectivity_2"
replace nombre = "Not selective program" if nombre=="Major_Selectivity_1"
replace nombre = "Selective program" if nombre=="Major_Selectivity_2"

drop 	if  regexm(nombre, "Selective") | regexm(nombre, "Not selective")
drop 	if  regexm(nombre, "P310")
drop 	if  regexm(nombre, "HS Higher Ed")
drop 	if  regexm(nombre, "Tercil") | regexm(nombre, "Stratum") 

replace p = 2 if regexm(nombre, "Male")
replace p = 3 if regexm(nombre, "Female")
replace p = 4 if regexm(nombre, "Not ethnic")
replace p = 5 if regexm(nombre, "Ethnic")
replace p = 6 if regexm(nombre, "First-gen college")
replace p = 7 if regexm(nombre, "Not first-gen college")
replace p = 8 if regexm(nombre, "Rural")
replace p = 9 if regexm(nombre, "Urban")
replace p = 10 if nombre == "HS SABER 11 Q1" 
replace p = 11 if nombre == "HS SABER 11 Q2"
replace p = 12 if nombre == "HS SABER 11 Q3" 
replace p = 13 if nombre == "HS SABER 11 Q4" 
replace p = 14 if nombre == "HS HQ college <P50" 
replace p = 15 if nombre == "HS HQ college >P50"
drop 	if p>15
labmask p, values(nombre)

twoway 	(rcap cs ci p if r=="11", lcolor(black)) ///
		(scatter cc p if r=="11", ysc(r(0)) mcolor(black) msymbol(O)) ///
		,ytitle(RD-IV Coefficient and 95% CI) yline(0) xlab(1(1)15,valuelabel angle(45)) ///
		legend(off) xsize(8) xtitle("") ylab(-.2(.2)1) ///
		caption(Sample restricted to SISBEN-eligible individuals, size(small)) graphregion(color(white)) 
		graph export "${outcomes}/FigC4a.pdf", replace width(20)

twoway 	(rcap cs ci p if r=="si", lcolor(black)) ///
		(scatter cc p if r=="si", ysc(r(0)) mcolor(black) msymbol(O)) ///
		,ytitle(RD-IV Coefficient and 95% CI) yline(0) xlab(1(1)15,valuelabel angle(45))  ///
		legend(off) xsize(8) xtitle("")  ylab(-1(.5)2.5)  ///
		caption(Sample restricted to SABER 11-eligible individuals, size(small)) graphregion(color(white)) 
		graph export "${outcomes}/FigC4b.pdf", replace width(20)
		
********************************************************************************
*** HETEROGENEITY: FIELD OF STUDY
********************************************************************************

use 	"${Intermediate}/heterogeneities_graduation_exam_pro.dta", clear

replace nombre = "cc" if _n==1
replace nombre = "cr" if _n==2
replace nombre = "sec" if _n==3
replace nombre = "ser" if _n==4
replace nombre = "cs" if _n==5
replace nombre = "ci" if _n==6
replace nombre = "cc0" if _n==7
replace nombre = "cr0" if _n==8
replace nombre = "pc" if _n==9
replace nombre = "pr" if _n==10
replace nombre = "bwc" if _n==11
replace nombre = "bwr" if _n==12
replace nombre = "nec0" if _n==13
replace nombre = "nec" if _n==14
replace nombre = "ner0" if _n==15
replace nombre = "ner" if _n==16
replace nombre = "n" if _n==17

xpose, clear varname

rename v1 cc
rename v2 cr
rename v3 sec
rename v4 ser
rename v5 cs
rename v6 ci
rename v7 cc0
rename v8 cr0
rename v9 pc
rename v10 pr
rename v11 bwc
rename v12 bwr
rename v13 nec0
rename v14 nec
rename v15 ner0
rename v16 ner
rename v17 n
rename _varname nombre
drop if _n==1

order nombre

gen 	o = "baseline"

foreach o in stem stem_plus arts social other {
	preserve
		use 	"${data}/heterogeneities_graduation_exam_pro_`o'.dta", clear
		replace nombre = "cc" if _n==1
		replace nombre = "cr" if _n==2
		replace nombre = "sec" if _n==3
		replace nombre = "ser" if _n==4
		replace nombre = "cs" if _n==5
		replace nombre = "ci" if _n==6
		replace nombre = "cc0" if _n==7
		replace nombre = "cr0" if _n==8
		replace nombre = "pc" if _n==9
		replace nombre = "pr" if _n==10
		replace nombre = "bwc" if _n==11
		replace nombre = "bwr" if _n==12
		replace nombre = "nec0" if _n==13
		replace nombre = "nec" if _n==14
		replace nombre = "ner0" if _n==15
		replace nombre = "ner" if _n==16
		replace nombre = "n" if _n==17

		xpose, clear varname

		rename v1 cc
		rename v2 cr
		rename v3 sec
		rename v4 ser
		rename v5 cs
		rename v6 ci
		rename v7 cc0
		rename v8 cr0
		rename v9 pc
		rename v10 pr
		rename v11 bwc
		rename v12 bwr
		rename v13 nec0
		rename v14 nec
		rename v15 ner0
		rename v16 ner
		rename v17 n
		rename _varname nombre
		drop if _n==1

		order nombre

		gen 	o = "`o'"
		tempfile `o'
		save "``o''"
	restore

	append using "``o''"
} 

keep 	if regexm(nombre,"male")	

* 		Plotting conventional RD coefficient and conventional confidence intervals
drop 	cs ci
gen 	cs = cc+(1.96*sec)
gen 	ci = cc-(1.96*sec)

gen 	r = "11" if regexm(nombre,"_saber11")
replace r = "si" if regexm(nombre,"_sisben")

gen 	p = _n
replace nombre = subinstr(nombre, "_saber11", "",.)
replace nombre = subinstr(nombre, "_sisben", "",.)
replace nombre = strproper(nombre)
labmask p, values(nombre)

gen 	x = .
replace x = 1 if o == "baseline"
replace x = 3 if o == "stem"
replace x = 5 if o == "stem_plus"
replace x = 7 if o == "arts"
replace x = 9 if o == "social"
replace x = 11 if o == "other"
replace x = x + 1 if nombre =="Female"

labmask x, values(nombre)

gen 	x2 = ""
replace x2 = "Any" if o =="baseline"
replace x2 = "STEM" if o =="stem"
replace x2 = "STEM Plus" if o =="stem_plus"
replace x2 = "Arts"  if o =="arts"
replace x2 = "SSH"  if o =="social"
replace x2 = "Other" if o =="other"

gen 	xx = x
labmask xx, values(x2)
replace xx = xx +.5

twoway 	(rcap cs ci x if r=="11", lcolor(black)) ///
		(scatter cc x if r=="11", ysc(r(0)) mcolor(black) msymbol(O)) ///
		(scatter cc xx if r=="11", xaxis(2) color(none)) ///
		,ytitle(RD-IV Coefficient and 95% CI) yline(0)  ///
		legend(off) xsize(8) ylab(-.05(.05).35) ///
		xtitle("Bachelor's Degree Earned by Field of Study", axis(2)) ///
		xtitle("", axis(1)) graphregion(color(white)) ///
		caption(Sample restricted to SISBEN-eligible individuals, size(small)) ///
		xlab(1(1)12, valuelabel axis(1)) xlab(, valuelabel labsize(normal) axis(2)) 
		graph export "${outcomes}/FigC5a.pdf", replace width(20)

twoway 	(rcap cs ci x if r=="si", lcolor(black)) ///
		(scatter cc x if r=="si", ysc(r(0)) mcolor(black) msymbol(O)) ///
		(scatter cc xx if r=="si", xaxis(2) color(none)) ///
		,ytitle(RD-IV Coefficient and 95% CI) yline(0)  ///
		legend(off) xsize(8) ylab(-.05(.05).35) ///
		xtitle("Bachelor's Degree Earned by Field of Study", axis(2)) ///
		xtitle("", axis(1)) graphregion(color(white)) ///
		caption(Sample restricted to SABER 11-eligible individuals, size(small)) ///
		xlab(1(1)12, valuelabel axis(1)) xlab(, valuelabel labsize(normal) axis(2)) 
		graph export "${outcomes}/FigC5b.pdf", replace width(20)	
